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Abstract 

We construct an approximate solution of field equations in the Abelian Higgs 
model which describes motion of a curved vortex. The solution is found to 
the first order in the inverse mass of the Higgs field with the help of the 
Hilbert-Chapman-Enskog method. Consistency conditions for the approxi- 
mate solution are obtained with the help of a classical Ward identity. We 
find that the Higgs field of the curved vortex of the topological charge n > 2 
in general does not have single n-th order zero. There are two zeros: one is 
of the (n-l)-th order and it follows a Nambu-Goto type trajectory, the other 
one is of the first order and its trajectory in general is not of the Nambu-Goto 
type. For \n\ — 1 the single zero in general does not lie on Nambu-Goto type 
trajectory. 
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1 Introduction 



In the recent paper |T|] we have shown how to construct approximate do- 
main wall solution with the help of Hilbert-Chapman-Enskog method H for 
solving equations with singular perturbations. The solution has the form 
of a perturbative series in positive powers of the width of the domain wall. 
Crucial role is played by certain consistency conditions for the perturbative 
expansion. Their appearence is related to existence of zero modes. In order 
to satisfy the consistency conditions we have had to allow for the possibility 
that the so called co-moving coordinate system is tied to an auxilliary mem- 
brane instead to zeros of the scalar field. Our approach differs essentially 
from earlier attempts to use the expansion in the width for domain walls. 

In the present paper we would like to apply this same method to obtain a 
solution in the Abelian Higgs model describing a generic non-selfinteracting 
vortex. By the non-selfinteracting vortex we mean a vortex such that no its 
parts which are distant along the vortex are spatially close to each other. 
This excludes, e.g., spikes and selfintersections. It turns out that the method 
works quite well also in the case of vortices. The calculations we report in 
the present paper are parallel to the ones in [|IJ. Therefore, we do not repeat 
here general motivation and comments about the course of calculations. In 
[p]] one can also find a more extensive list of references including also papers 
on dynamics of vortices. 

The vortex in the Abelian Higgs model is the classic example of a vortex. 
It is called local because its fields are different from their vacuum values 
only in a close vicinity of a single line - far from this line the fields are 
exponentially close to the vacuum values. Thus, there are no long range forces 
acting on pieces of the vortex. For this reason the dynamics of such a vortex 
is probably simpler than dynamics of a global vortex with its long range 
Goldstone field. This is the reason we have chosen to apply the new method 
first to the local vortex. Unfortunately, the Abelian Higgs model is not simple 
on the mathematical side. It involves six real fields (ReQ, Jm$, A^) which 
obey a non-linear set of equations. Moreover, except for the special case of 
Bogomol'nyi limit there are present two length scales, given by the inverses of 
masses of the scalar and vector fields. Even in the simplest case of a straight- 
linear, static vortex the exact analytical form of the corresponding solution 
is not known. Therefore, one should be prepared for rather cumbersome 
calculations. 
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The plan of our paper is the following. In the next Section we introduce 
the co-moving coordinate system and we write the field equations in the new 
coordinates. Most of this material is not new. We quote it here for conve- 
nience of the reader and in order to fix our notation. Moreover, we would 
like to emphasize the fact that the co-moving coordinate system apriori is 
not tied to position of zeros of the Higgs field. The approximate position 
of the zeros are determined at the end, from the constructed approximate 
solution of the field equations. In Section 3 we find the zeroth order term in 
the expansion in the width of the vortex. By the width we mean the inverse 
mass of the Higgs field. We also check the consistency conditions appropri- 
ate for the zeroth order and characteristic for the Hilbert-Chapman-Enskog 
method. Section 4 is devoted to the first order solution and the relevant 
consistency conditions. Section 5 contains final remarks. In Appendix A we 
derive a classical Ward identity which helps to calculate troublesome inte- 
grals encountered in Section 4. In Appendix B we give a list of constants 
appearing in Section 4. 

2 Equations for the scalar and vector fields 
in the co-moving coordinates 

We take the Lagrangian of the Abelian Higgs model in the following form 

1 A M 2 

L = DWDp* - -F^F, V - -(|<t>| 2 - — ) 2 , (1) 

where 

The Minkowski space-time metric is chosen as (r)^ u ) = diag(l, — 1, — 1, — 1). 
The corresponding field equations are 

A M 2 
ZW^ + -$(|$| 2 -— ) = , (2) 

& i F llv = iq(&D v &-QD v &). (3) 
The vacuum fields are characterised by the following conditions 

|$|=$ = -^=, D^ = 0. (4) 
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The masses of the Higgs and the vector particle are equal to m 2 H = M 2 , m\ = 
^-rnjj, respectively. We consider the case ^f- < 1 which corresponds to 
superconductors of the second kind. 

In the case of single vortex considered in this paper, the fields $, at 
each instant of time have the vacuum values in the whole space R 3 except for 
a vicinity of a line C where the energy density is significantly different from 
zero. We do not assume that the Higgs field vanishes precisely on this line. 
The world-sheet of the line C in Minkowski space-time we shall denote by 
S. The line C we will call the string, in anticipation of the result of the next 
Section, where we find that in the framework of 1/M expansion C has to obey 
an equation identical in its form with Nambu-Goto equation for a relativistic 
string. The string is used to define the coordinate system co-moving with 
the vortex. To this end, we parametrize S with two parameters r, a 

[y"(r,<7)]eE. 

The parameter r is by assumption time-like, i.e. Y^Y^ T > 0, while a is the 
space-like one, i.e. Y^Y^ C < 0. In the following we use a more compact 
notation: (r, a) = (u a ), a = 0, 3, with a = (a = 3) corresponding to r (cr). 
The co-moving coordinates are defined in a vicinity of the string world-sheet 
E by the following formula 

x » = Y fl (u a ) +C <(w a ), (5) 

where x M are Cartesian, laboratory frame coordinates in Minkowski space- 
time; i = 1, 2; ni(u a ) are two orthonormal space-like four-vectors, orthogonal 
to E at the point [Y M (ii a )] in the covariant sense, i.e. 

Y£(u) n ifl (v) = 0, nt{u)n klJb {u) = -5 ik . (6) 

The conditions (6) do not fix the four-vectors rii uniquely - there is still a 
freedom of local (i.e. (w a )-dependent) SO (2) rotations acting on the indices 
i, k. The coordinates £ x ,£ 2 vanish on the string. We shall denote the full 
set of co-moving coordinates by (C°0 : (C°0 — ( ua ,€), a=0,3,l,2. The 
co-moving coordinate system has been used in numerous papers, see, e.g., 
[@- [|, |> Uli with the difference that in those papers it was attached to a line 
on which the Higgs field was assumed to vanish. 
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The next step is to write Eqs.(2),(3) in the co- moving coordinates. It is 
convenient to introduce several mathematical objects: 

- the metric g ab on £, 

9ab — Y^Y^bi 

- the extrinsic curvature coefficients, 

Kab — ~ n i^,ab', 

- the torsion coefficients 

where e tk is the 2-dimensional totally antisymmetric symbol (e 12 = 1). The 
torsion and extrinsic curvature coefficients depend on the choice of the four- 
vectors rii. In general it is not possible to adjust them in such a manner that 
the torsion vanishes. This can be seen from Ricci equation 

d a uo b - d b uj a = e ij K a c K J bc , 

where K l a c = K ab g bc . Vanishing of the torsion would imply that the extrinsic 
curvatures at a given point could not be arbitrary. 

The Minkowski space-time metric tensor in the new coordinates has the 
following form 

(7) 
where 

c ah = 9ab + 2Ki b c + K^ee - ee^aco b , 

G a i = Gi a = e tk £ h uj a , 
Gik = —$ik- 
It follows from these formulae that 

where G = det[G a p], g = det[g a b], and 

h(u a , e) = i + kx + \{K a Ki b - K:,»Ki")ee. (8) 



[Gap] 



G„ 



Gi a Gi 



G a k 

ik 
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We shall also need components of the inverse metric tensor They are 

given by the following formulae: 

G ab = h- 2 [g ab (l + C Kl c f - 2CK iab {l + i k K k c c ) + i' i' K'," K jl " . (9) 

G m = e lk i k u b G ab , (10) 
G ik = _ S ik + (^fc^^i _ £ i £ k ) U j a uj b G ab . (11) 

We shall write the Higgs and Maxwell equations (2), (3) in the co-moving 
coordinates using rescaled variables: 

C = jjrS\ $K,f) = $ O 0K,^). 

Also the vector field is transformed to the new coordinates. It has the 
following components 

(A a ) = (A a , Ai), 

where a = 0, 3, i — 1,2, and 

a/2 

A a = e£(u,0A li , Ai — —rtf(u)A^. 

Here 

Notice that the formula defining the components Ai contains the factor 
v^2/M, so that Ai are dimensionless, similarly as s l and 0. We do not rescale 
the coordinates u a and the components A a . The variables £ l present in G a/3 
and \/—G are also rescaled by the same factor. We shall use the following 
notation for the U(l) gauge covariant derivatives 

d d 
D a = — + tqA a , Di = — + iqAi. 

Using the rescaled co-moving coordinates and the rescaled fields we can 
write Eq.(2) in the following form 

AA0 + 0-|0| 2 = -^A0 (12) 
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(D a + f?s l u a Di) ^hG ab (D b <j) + e kr s r u b D k <t>) 



ab i 



M 2 ^gh 



The function h is given by formula (8); denotes d/ds\ 

The Maxwell equation (3) in the co-moving coordinates has the form 



(13) 



where 



dA^ dA 



Pi 



d(? do ' 

It is convenient to split Eq.(13) into two equations and to pass to the rescaled 
variables. In the case 7 = k — 1, 2 Eq.(13) gives 



+ 



2A 



M 2 ^gh 



Xd t F lk + 2iq(<j>*D k (l> - 0/^0*) = ~^F lk 



{d a + e l s l u a di 



(14) 



where 



Fi k — diA k — d k Ai, F bk 
For 7 = c = 0, 3 we obtain from (13) 



-ghG ab (F bk + e sr s r uj b F sk ) 
d b A k - d k A b . 



XdiF ic + 2iq{(t>*D c (t> - 0D C 0*) = ~^F ic 



(15) 



2A 



M 2 ^/=g~h 



(d a + e tk s k u a d t ) V=ghG ab (F bc + e tl s l u b F tc ) 



Equations (12), (14), (15) have the form convenient for constructing the 
perturbative expansion in powers of 1/M. 

The reason for the splitting of the Maxwell equations into the two groups 
is that the components A± play different dynamical role in the motion of 
the vortex than A c . The components A^ determined basically from Eq.(14), 
do not vanish for any vortex and they are present even in the zeroth order 
approximation. The components A c are determined from Eq.(15). They 
vanish in the zeroth and first order approximations for all vortices, and for 
some vortices to all orders (e.g. for a straight-linear vortex). Moreover, we 
shall see that in the 1/M expansion we first find the n-th order contributions 
to (j> and Ai and the n-th order contribution to A c is calculated afterwards. 
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The co-moving coordinate system in general does not cover the whole 
space-time. We assume that the range of its validity includes the physically 
most interesting region - outside of it the fields are exponentially close to 
the vacuum. The range of validity of the co-moving coordinates can be 
determined from the condition h(u a ,^ 1 ) > 0. We will not discuss here this 
purely mathematical point. 



3 Expansion in the width — the zeroth order 
solution and zero modes 

The present Section is devoted to a preliminary discussion of the expansion 
in the width for the vortex solution constructed with the help of Hilbert- 
Chapman-Enskog method. The zeroth order solution we obtain coincides 
with the one known from earlier applications of the expansion in the width 
to vortices, 0. In the next Section we calculate the first order contribution. 

We seek approximate solutions of Eqs.(12),(14),(15) in the form of power 
series in 1/M: 

i ... i 

(16) 

(17) 
(18) 

We also assume that these fields do not have any component oscillating in r 
with a frequency proportial to M. This assumption ensures that when act- 
ing with the derivative d T on the n-th order contribution we obtain still the 
n-th order contribution. To illustrate the point: the r-derivative of the n-th 
order function M~ n cos Mr belongs to the (n-l)-th order. This assumption 
has a physical content - it eliminates a whole class of solutions, in particular 
the ones describing radiation of the Higgs and vector particles. In order to 
include such solutions one would have to construct an extended perturbative 
scheme, e.g., in analogy to considerations presented in f?|, where such oscil- 
lating component was taken into account in the case of domain walls in the 
framework of a polynomial approximation. 

In the perturbative calculations we would like to preserve the local U(l) 
gauge invariance of the Abelian Higgs model. To this end, we use the trick 



= 




+ -</> (1) - 
M 


^ 2) + - 


M = 




+ 1 A {1) 




A a = 


4 0) 


+ M A « ] 


+ M^ 2) + 
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borrowed from background field techniques in field theory: we assume that 
the gauge transformations act only on , while with n > are gauge 
invariant. Then, the total vector field A a has the right transformation law. 
As for the Higgs field, each contribution <f>( n ) is multiplied by the same (C )- 
dependent phase factor. Then, in the gauge covariant derivatives it is suffi- 
cient to use the zeroth order vector field i.e. 

= + DP* = ( A + 

As explained in detail in paper the expansion in 1/M is an example of 
singular perturbations. In order to calculate the perturbative contributions 
correctly we use the Hilbert-Chapman-Enskog method. In this method, apart 
from obtaining the solution in a given order it is essential to look also at the 
equations in the higher orders - from them one obtains consistency conditions 
for the solution in the given order. The consistency condition for the zeroth 
order solution follows from the equations in the first order. The equations 
in the zeroth and first orders are obtained by inserting the expansions (16- 
18) into Eqs.(12),(14),(15) and collecting the terms of the order (1/M)° and 
(1/M) 1 . 

In the zeroth order, Eqs.(12),(14),(15) are reduced to the following equa- 
tions: 

D (°>(£>(°y0)) + (O) _ | (O)|2 (O) = 0> (1Q) 

= -2.g[0*(°)(^ o) (O) ) - ^(4V 0) )1, (20) 
XdiF® = -2iq[^\Df^^) - 0(°)p(°)0(o))*]. (21) 

Equations (19), (20) coincide with equations for the static Abrikosov-Nielsen- 
Olesen (A-N-O) vortex §, [§. 

Let us quote here basic facts about the A-N-0 vortex. We shall use them 
in the next Section. The axially symmetric Ansatz for the A-N-0 vortex has 
the form 

k 

0(0) = e M> F , 8 ) AO) = t ik s _^ _ H , 8 ss (22) 

s qs 

where s = J (s 1 ) 2 + (s 2 ) 2 , 9 is the azimuthal angle in the (s 1 , s 2 )-plane, and 
n is the topological charge of the vortex. This Ansatz reduces Eqs.(19),(20) 
to the following equations 

F" + -F' - q 2 H 2 F + F - F 3 = 0, (23) 
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B'(s) = -j-F 2 H, (24) 

where 

B = H'+^ 

s 

and ' denotes the derivative d/ds. These equations are supplemented by the 
boundary conditions for the functions F, H: 

F(0) = 0, F(oo) = 1, (25) 

71 

limsH(s) = - H{oo) = 0. (26) 

We will call the functions F,H obeying Eqs.(23), (24) with the boundary 
conditions (25), (26) the A-N-0 functions, and denote them by F (s), H (s). 
They exponentially vanish for large s. The magnetic field of the vortex is 
equal to 

FM = e ^ + H ')=e ik B (s). (27) 
Explicit, exact analytical form of the A-N-0 functions is not known. Never- 



theless, they are very well investigated, see, e.g. |L0 . 

For 0^,t4j-°' ) given by the A-N-0 solution, equation (21) is reduced to 

A4 0) - ^(s)A^ = 0, (28) 

where A = didi is the 2-dimensional Laplacian. The operator on the l.h.s. 
of this equation does not have non-zero regular, localised solutions. This can 
be easily seen by interpreting the operator —A + -%-Fq as a Schroedinger 

operator with the potential ^Fq. Because this potential is non-negative, 
there are no bound states with zero eigenvalues. Therefore, 

= 0. (29) 

Notice that Eq.(28) coincides with Eq.(24) for H(s). The fact that the solu- 
tion Hq(s) of Eq.(24) is non-zero is not in contradiction with the conclusion 
(29) because H is singular at s = (see (26)), while A^> is required to be 
regular for all s l and exponentially vanishing for large s. These regularity 
conditions follow from the requirement that the fields F a p should be regular 
and exponentially localised. 
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Notice that by adopting the concrete form of the zeroth order solution 
we implicitely fix the gauge for two local gauge symmetries. The first local 
gauge group is formed by the (x M )-dependent "electromagnetic" U(l) trans- 
formations. The other one consists of (w a )-dependent SO (2) rotations in the 

£ 2 )-plane, equivalent to SO(2) rotations of the basis (rii,n 2 ) in the plane 
orthogonal to the world-sheet S of the string. The third local gauge symme- 
try, namely invariance with respect to reparametrisations of the world-sheet 
of the string remains explicitely preserved. In the present paper we prefer 
to work with the concrete gauge fixing for the former two gauge symmetries. 
Actually, one can provide a formulation of the Abelian Higgs model in the 
co-moving coordinates which is explicitely invariant under the three local 



gauge symmetries, [III 



Now we will turn to the consistency conditions for the zeroth order so- 
lution. To this end we have to derive equations in the first order. It is 
convenient to introduce a compact notation for the fields and the field equa- 
tions. Thus, 

* = {^t \ ■ (30) 

The compact notation for the left hand side of the set of equations formed 
by Eq.(12), its complex conjugate counterpart, and Eq.(14) is J 7 ^). Then, 
the zeroth order equations (19), (20) can be written as 

= 0, (31) 

/ 



0*(o) 

V4 0) 



(32) 



In this notation we regard and 0* as two independent fields obeying Eq.(12) 
and its complex conjugate countertpart, correspondingly. This is equivalent 
to considering the real and imaginary parts of the field <fi and Eq.(12). Notice 
that we have left aside Eq.(15) - it will be solved later. This is possible 
because A a is not present on the l.h.s. of Eqs.(12),(14). 
Let us now consider a small correction 



5V = 5<p* (33) 
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to the zeroth order solution Because of Eq.(31) we have 

F(yW + 6*)=L6* + 0((6*)). (34) 

Here L is a linear operator obtained by expanding the l.h.s. of Eqs.(12),(14) 
around (f)^°\ Af\ Including also 0* and the complex conjugate of Eq.(12) we 
find that 



L = 



D f) jD f ) + l-2|0(°)| 2 

_(0*(O))2 



d *(o) d: (o) + l 



_(0(O))2 

-2|0(°)| 2 
-2iq[^D^ - (D<V 0) )] 

2ig( J D{ O) (O) ) +iq<P {0) di 
-2iq(Dl 0) <f> {0) y-iq<f>* {0) di 



.(35) 



The equation for the first order correction to \E'^ - ) , i.e. for 



v4" 



(36) 



easily follows from formula (34) and Eqs.(12),(14). It has the following form 



/ D 



(0) 



= -y/2 



(0) 
ik 



Kl 



(37) 



where the r.h.s. has been obtained by expanding the r.h.s. of Eqs.(12), (14). 
Equation (37) will be solved in the next Section. 

From Eq.(37) we obtain the consistency conditions for the zeroth order 
solution. The crucial observation is that there are left eigenvectors of operator 
L with the corresponding eigenvalue equal to zero, i.e. normalizable, non-zero 

■fyl.z.m. guch that 

' ds 1 ^^ 12 '" 1 -)*^* = (38) 



for all \& vanishing sufficiently quickly for large s. In condition (38) T denotes 
the matrix transposition, and the abbreviation l.z.m. stands for "left zero 
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mode" . It is easy to check that the operator L is Hermitean with respect to 
the scalar product given by 



< ^1*2 >= J ds l ds 2 m* T m, 



(39) 



in the space of Ws which have the second component equal to the complex 
conjugation of the first one. (Therefore, it is a Hilbert space over real numbers 
and not over complex numbers.) Actually the Higgs self-coupling constant 
A has been placed in front of the F a p field in Eqs.(14),(15) just in order to 
obtain L Hermitean with respect to the simple scalar product (39). In the 
case of Hermitean operator the l.z.m.'s coincide with right zero modes, which 
are defined as normalizable (non-zero) solutions of the equation 



= 0. 



(40) 



The zero modes can be easily found with the help of formula (34) by 
using a particular 5^ corresponding to symmetries of Eqs.(19),(20). These 
equations are regarded here as equations for functions on the plane (s 1 , s 2 ). 
Because of invariance of Eqs.(19),(20) with respect to translations in the 
(s 1 , s 2 )-plane, one can generate from the solution infinitely many other 
solutions by applying the translations. We shall denote these new solutions 
by T *(°). Thus, 

/ 0(°)(s* + a 1 ) 



M (o) 



0*(O) ( 

v 4 0) ( 



where a = (a 1 ) is a constant vector defining the translation. It is clear that 

F(T a ¥ 0) ) = 0. 
Taking infinitesimal a, we obtain from formula (34) that 

L% = 0, 



where 



d(T a *W) 



d<p* (0) 

ds l 



dA 



(») 



V ds l ) 
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These can not be used as the zero modes because they are not normal- 
izable. The correct zero modes are obtained by combining the translation 
with a gauge transformation in the spirit of improved variations of Jackiw 
and Manton, 11121. Instead of T n ^^ we take 



f *(o) 



s l + a l 



(41) 



where x{ a 'i s *) = qcL l A[°\s l ). It is clear that also T a ^^ is a solution of 
Eqs.(19),(20). As the zero modes we take 



_ _ d(f a *<®) , 



I A ( V 0) \ 
(£>f°)0(o))* 



V 



,(0) 
Ik 



(42) 



/ 



where I = 1,2. They are normalizable and orthogonal - one can easily check 
that 

< e fc |0 z >= S kl J ds^s^Bl + Fo 2 + q 2 H 2 X), 

where F , H , and B are the A-N-0 functions. 

There also exists a solution of Eq.(40) related to the U(l) gauge invari- 
ance. It has the form 



N(0) 



V 







Another possible solution, corresponding to rotations around the origin in 
the (s 1 , s 2 )-plane is in fact identical with \l/ 5 because in the case of the A-N- 
O vortex the rotations are equivalent to global U(l) gauge transformations. 
Notice however that is not normalizable with respect to the scalar product 
(39). Therefore, it can not be accepted as the zero mode. 

The consistency conditions follow from Eq. (37) by integrating it with the 
zero modes 0;, like in formula (38). The r.h.s. of Eq.(37) in general is not 
a linear combination of the zero modes because of the factor A in front of 
Nevertheless, the integrations are simple and we obtain the following 



? (o) 

Ik ■ 



conditions 



K 



la 



(43) 
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where I = 1, 2 and a = 0, 3. 

The conditions (43) are equivalent to Nambu-Goto equation for a rela- 
tivistic string whose world-sheet is given by the surface £ in Minkowski 
space-time. 

To summarize, the zeroth order solution is given by formulae (22) with 
F(s), H(s) equal to the A-N-0 functions F (s), H (s). The components 
vanish. The variables s = M£* are related to the Cartesian coordinates 
by formula (5). For the functions Y >1 (u a ) in that formula we can take any 
solution of Nambu-Goto equation (43). In the zeroth order approximation 
the string coincides with the zeros of the Higgs field. 



4 The first order corrections 

The first order corrections to <p an d obey Eq.(37) with vanishing the 
r.h.s., due to the consistency conditions (43). This equation is satisfied by 
linear combinations of the zero modes and with coefficients which do not 
dependent on , i.e. 

*W = C'K)0 z (s l ) + C 3 {u a )V g {s l ), (44) 

where I = 1,2. Now we have included ty g because in Eq.(37) there are no 
integrals. The functions C l ,C 3 are real, otherwise the second component of 
of would not be equal to the complex conjugation of the first one. In 
the following we shall show that the functions C , C 3 are not arbitrary. 

We also have to solve the first order equation obtained from Eq.(15) for 
Ab. With the Nambu-Goto equations (43) and the zeroth order results taken 
into account, it has the following form 

A4 1} - ^(s)A^ = 0. (45) 

We already know from the discussion following Eq.(28) that the only regular, 
localised solution of Eq.(45) is the trivial one, 

AP = 0. (46) 

By the analogy with the case of domain walls [|l|] we expect consistency 
conditions for the first order solution coming from the second and third order 



14 



equations. These conditions have the form of equations which must be obeyed 
by the functions C l ,C 3 . Therefore, we have to consider the perturbative 
equations in the second and third orders. As for still higher orders, we 
expect that the corresponding consistency conditions can be "saturated" by 
new arbitrary functions, introduced in each order as the general solution of 
the homogeneous equation L\&( n ) = 0. All such solutions have the form (44) 
with the functions C l , C 3 replaced by a new set of functions in each order. 
They have to be added to a solution of the full perturbative equation in 
the given order to obtain the general solution in the n-th order. The full 
perturbative equations in the second and third orders are given below, see 
Eqs.(47),(52). 

The equation in the second order has the following form 



( / (2) 

J*(2) 

\ 4 2) 



(47) 



where 



/( 2 ) = 20(°)|0( 1 )| 2 + (0( 1 )) 2 0*(°) + g 2 Aj 1) 4 1) (O) 
-2iqA^ A ( °V (1) " iqidiA^)^ + 2K? K^sWf } (o) 



(48) 



(d a + e il s l u a Df ) )i < y/=g S f b ^ r s r u b D 



,(o) 



and 



a< 2 > = -2i q (^D^ W - 0«4V (1) ) 
+2g 2 4 1) (0* {1 ¥° ) + <f>m<f>*W) + 2\KfKl a s^F^ 
+ 



j±(d a + e u s l cu a d t )(^g ab e tr s r LU b F t f) 



(49) 



We also have the following equation for Af> (obtained by expanding Eq.(15)) 



(2)> 



a4 2) - T^o(^4 2) = d c (d t A\ 

_2M(0*(O),9 c 0(2) _ 0(0)^(2) + 0*(l)9 c 0(l) _ d c (/>*W). 



(50) 



In the formulae for f^ 2 \a^\ and on the r.h.s. of Eq.(50) the zeroth order 
results and the Nambu-Goto equation (43) have already been taken into 
account. 
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From Eq.(47) one can in principle determine the second order correction 
. In practice we do not hope to obtain an explicit formulae for - we do 
not even know the explicit form of the zeroth order functions (f)^ , which 
are present in the operator L. On the other hand, the functions 4>(°\ are 
known numerically, and the coefficients and the r.h.s. in Eq.(47) are regular. 
Therefore, we expect that it will be possible to investigate the solutions of 
this equation by numerical methods. This task we will leave for a separate 
work. In the present paper we would like to focus on the description of the 
perturbative scheme and on the first order correction 

As explained in Section 3, the consistency conditions are obtained by 
integrating the both sides of the perturbative equation, i.e. Eq.(47) in the 
present case, multiplied by the zero modes 6^: 

= / d*W[(D\ O) (O) )7 {2) + c.c. + F^a?] 

(c.c. stands for "complex conjugate"). Most contributions to the integral 
vanish because the integrands are odd functions of the variables s l and the 
integration range is from — oo to oo. The non-vanishing terms on the r.h.s. 
lead to the following condition 

= -d q e lk C l (u a )C 3 (u a ), 

where I, k = 1,2, and d is a non- vanishing constant given in the Appendix 
B. 

We shall see that in general C l are not equal to zero. Therefore, we put 

C 3 = 0. (51) 

Equation (50) does not yield any restrictions on the functions C l , C 3 
because the operator on the l.h.s. of it (the same as in Eq.(28)) does not 
have zero eigenvalues. 

Further consistency conditions follow from the equations in the third or- 
der. These conditions have the form of non-linear, intercoupled, inhomoge- 
neous wave equations for C l (u a ) regarded as fields on the world-sheet S of 
the Nambu-Goto string. Unfortunately, now the calculations become rather 
cumbersome - one could hardly expect simple formulae in the third order 
approximation in the model with six coupled fields. The equation for the 
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third order correction has the form 



( / (3) 

/*(3) 

I 4 3> 



(52) 



where 



/ (3) = A (3) + / : 



r(3) _(3) 



l {3) 



+ a 



(3) 
2fc • 



(3) 



a u! ( /i 3 "*? a 2fc^ ) denote the terms obtained by expanding the 



Here fl 

l.h.s.'s ( the r.h.s.'s ) of Eqs.(12), (14), respectively. We find that 



-2iqAf ] Df V (2) - 2iqAf\Df + ^4 (1) (O) ) (53) 



af fc ) = -2i g [0*( 2 )4 o V (1) - 0( 2 )(4 O) (1) )* 

+<j)*WD^<j) {2) ~ (1) (4°V (2) )*] (54) 
+4g 2 4 1) (0* (2) (O) + { V (O) ) + 4g 2 4 2) (0* (1 V (O) + (1) 0* {O) )- 

and 



+2ige l 's*e fcr s r w a ^ a6 ( J Df ) {O) )4 1) (55) 
f ^=(9 a + e l 's^ a A (0) )[v / =^^(^ {1) + e fcr s r ^( J Di° ) (1) + ig4V 0) ))] 



f ^(9 a + e <i s'w a ^)[V=^^(F l £ ) + e'Yw^?) (56) 
-2V2 v ^gKP ab to b s p e tr s r F^ ) ] . 
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In these formulae we have taken into account the Nambu-Goto equation as 
well as the fact that vanishes. 

In the third order we also have an equation for . We will not present 
it here because it does not lead to consistency conditions - the operator on 
the l.h.s. of it is the same as in Eqs.(28),(45). 

The next step is to integrate over s 1 ,^ 2 the r.h.s. of Eq.(52) multiplied 
by the zero modes. At first look this might seem a difficult task because the 
second order corrections <f)( 2 \ A^ are given only implicitely by Eq.(47). The 
other functions, including the zero modes, are explicitely expressed by the 
A-N-0 functions Fq(s),Hq(s). Luckily, it turns out that the 1/M expansion 
has the special property that the troublesome integrals can be expressed 
by integrals involving only the zeroth and first order corrections, which are 
explicitely given in terms of the A-N-0 functions. Derivation of this relation 
is based on the translational and U(l) gauge invariances of the model, and 
it strongly reminds derivations of Ward identities for Green functions in 
quantum field theory. In the quantum case, instead of the integration over 
s l ,s 2 there is a functional integral over fields. For this reason, we will call 
the identity presented below the classical Ward identity. 

All the terms with <f>( 2 \ <f>*( 2 \ present on the r.h.s. of Eq.(52) are 
obtained by expanding the l.h.s. 's of Eqs.(12),(14). They are collected as 
A (3) and 

aj^ • The crucial observation is that they follow from the expression 

v M ' 

by expanding it with respect to 1/M and keeping the term linear in 1/M 
(with 1/M dropped out). The operator L{^!^> + jj^^) is given by formula 
(35) with *(°) replaced by ^ (0) + The integrals we are discussing have 

the form 



8(1/ M) 



¥ 2 \ 

1/M=0 



where 1 — 1,2 and is given by formula (44) with C 3 = 0. In the Appendix 
A we prove that 

( D^Df^) V 

J, = / ds 1 ds 2 C r 



_d_F (0) 

ds r Ik 



L(¥^)¥ 2 \ (57) 
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The expression L(^/^)^/^ on the r.h.s. of the classical Ward identity (57) 
is eliminated with the help of Eq. (47) . In this manner X\ is given by integrals 
involving the zeroth order solution and the first order corrections only. 

Further steps in obtaining the consistency conditions are nothing more 
than laborious calculating of the integrals present in the initial form of these 
conditions which is 



0=li + j ds'ds^DfU^Tff + c.c. + F^a. 



(3)i 
2k J 



The result we have obtained has the form of the following equation for C l 
U^C 1 - 2cu a g ab e lk d b C k - u a u b g ab C l - -±=d a {^g ab u b )e lk C k 

- K lb K{ a C 3 + ^uj a g ab (e lk d h C k + uj b C l ) - ^(C k C k )C l = (58) 

d 1 yj-g di 



where 



□ (3) C*' = -^=d a {^g ab d b C l ). 



The constants di, i=l,...5, are listed in Appendix B. The fact that the coeffi- 
cients in Eq.(58) are correlated is due to the implicit local SO(2) invariance 
mentioned below formula (29). 

Let us summarize our results. The obtained to the first order in 1/M 
vortex solution has the following form ( in the co-moving coordinates ) 

= e in9 [F (s) + ^(u^jF^s) - iqe ikS ^H (s)F (s))}, (59) 

A = e ik j(^ - H (s)) - -L C k ( U y k B (s), (60) 

A b = 0, (61) 

where C k {u) obey Eq.(58) and B = H' + H /s. The problem how to fix 
initial data for equation (58) is discussed below. 
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From formulae (59), (60) one can find approximate position of the zeros 
of the Higgs and vector fields. Using the fact that for small s l [|K| 

where Cq is a constant, we find that the n-th order zero of at s 1 = s 2 = 
is split by the first order correction into (n-l)-th order zero at s 1 = s 2 = 
and a first order zero at 

4 « ~^C\ (62) 

For \n\ = 1 only the zero (62) is present. Thus, the (n-l)-th zero lies on the 
string and therefore follows a Nambu-Goto type trajectory. The other zero 
is shifted from the Nambu-Goto trajectory if the C l functions do not vanish. 

As for the zeros of the Ak field (k=l,2) , we find that the both components 
vanish at the point 

Sa « ~^C\ (63) 

which is in general different from the origin and s l H . This first order zero is 
not so important because its position can be shifted by the local U(l) gauge 
transformations dependent on s 1 ,^ 2 . 

The final problem we have to deal with is to fix the initial data for Eq.(58) 
from which we calculate the functions C l , and also for the Nambu-Goto 
equation (43) from which we find the functions Y^(u a ). The input is the 
initial data for the vortex fields in the laboratory frame coordinates. The 
string is only an auxilliary mathematical object. Therefore, we can choose 
its initial position and velocity in a convenient manner. 

For |n| = 1 the most natural choice is that at the initial r = To the 
position and velocity of the string coincide with the ones of the line of zeros 
of the Higgs field. This implies that at r = tq 

C l = 0, d T C l = 0. (64) 

The values of the C l functions at later time depend on the r.h.s. of Eq.(58). 
Only if it vanishes for all r then C l remain equal to zero for all r. Notice 
that for non- vanishing of r.h.s. of Eq.(58) it is necessary that the torsion 
coefficients u a do not vanish identically. 

For \n\ > 1 the initial data can be fixed again by looking at the zeros 
of the initial Higgs field. In general, it is possible that the zeros are split 
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already at the initial time and that they have nonvanishing relative velocity. 
For \n\ > 2 the split zeros are of different order - the initial position and 
velocity of the string are determined from the (n-l)-th order zero, while the 
initial values of the functions C l ,d T C l are fixed by the initial position and 
velocity of the line on which lie the first order zeros. For \n\ = 2 we can 
take any of the two zeros to determine the initial data for the Nambu-Goto 
string. Then, the other one will give the initial data for C l . The two choices 
are equivalent because Eq.(58) is symmetric with respect to the change of 
sign of the functions C l . 

In all cases the functions C l (u a ) give deviations of the trajectory of the 
first order zero of the Higgs field from the Nambu-Goto trajectory. 

5 Remarks 

In the present paper we have adapted the Hilbert-Chapman-Enskog method 
to the problem of evolution of the vortex. The method, previously applied 
to domain walls seems to work quite well also in the case of vortices, 
eventhough in the present case the calculations are more complicated. 

There are important differences between the aproach proposed in our 
paper and the earlier applications of the expansion in the width to dynam- 
ics of vortices in Minkowski space-time, 0, |], |5|, |J. Apart from using the 
Hilbert-Chapman-Enskog method in which the consistency conditions play 
the prominent role, also new is the fact that the co-moving coordinate system 
does not have to be tied to the zeros of the Higgs field. 

Already the first order correction reveals interesting phenomena in 
the vortex dynamics, like the splitting of the zeros of the Higgs field, or the 
possibility that for the \n\ = 1 vortex the zeros do not follow the Nambu- 
Goto type trajectory. It would be also interesting to know the second order 
correction, but it can be calculated probably only by means of numerical 
methods. 

Detailed discussion of properties of the obtained vortex solution, as well 
as concrete examples of evolution of the vortex obtained by solving Eqs.(43), 
(58) we will present in a separate paper, because the present one already 
seems to be rather voluminous. 

Let us end the present paper with a general remark that in our approach to 
the vortex dynamics there is a quite interesting interplay of the 4-dimensional 
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and 2- dimensional field theoretical models. First, the zeroth order solution 
as well as the zero modes are determined from Eqs.(19),(20) which can be 
regarded as equations of the 2-dimensional Euclidean Abelian Higgs model. 
Second, the functions C l (u) can be regarded as fields on the 2-dimensional 
Nambu-Goto manifold S. This non-trivial background manifold is to be 
determined from the Nambu-Goto equation (43). We expect that in higher 
orders new sets of functions on E will also appear. In this manner, the 
original description of the vortex by the 4-dimensional fields $, would be 
transformed into another description in terms of the Nambu-Goto string and 
the fields on its world-sheet. 



APPENDIX A. The derivation of the classical Ward identity (57) 



We start from the identity 

o = J ds'ds 2 (f a e fc ) t L(f a ¥^) f a ¥ 2 \ (Ai) 

valid for any a independent of s^s 2 . The operator T a is given by formula 
(41). The identity (Al) is just Eq.(38) written in terms of the translated and 
gauge transformed fields. It is valid for any \p( 2 ) with the proviso that the 
integral is convergent, but here is of course the second order correction. 
We take the particular a equal to 

, CHu) 



M 



where C l are the functions appearing in the first order solution (44). Let us 
recall that C 3 = 0. Next, we expand in 1/M: 

f $(o) ~ ^(o) , J_^(i) 



A (0) 













D *(o) 













d 

w . 
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and 



7^(2) ^ ^(2) + i^( 2 ) 



(the precise form of is not needed). Substituting these formulae on the 
r.h.s. of identity (Al) we obtain up the first order in 1/M 

o = J ds l ds 2 [(e fc )t ^( 2 ) + -^(e fc )t a* < 2 > 



M J 



ds 2 C l 



Df ] 



*(0) 



o o 



o 



V 



L(#(°)) *( 2 ) 



1 1 

M fc .1/2' 



The first integral on the r.h.s. vanishes because 0^ is the zero mode, and we 
obtain the identity (57). 



APPENDIX B. The constants appearing in Section 4 

In order to obtain the formulae listed below we have used Eqs.(23),(24), 
the boundary conditions (25), (26) and we have applied multiple integration 
by parts. 

d = 2rc dssB (F 2 - F 4 + F 2 + q 2 H 2 F 2 ), 
Jo 



d 1 = / dss(XB 2 + F 2 + q HqFq), 
Jo 



9 r°° 
d2 = -A / dssBl, 
2 Jo 



ds = \ Jo°° + (§- 4 -£)(qH F y + (1 - %){qFlH F 

„4 



+ (2qB ) 2 (F^ 2 + q 2 H 2 F 2 ) + 2 -^H 2 F^l - F 2 ) + ^H^F^H, + B )}, 
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di = V2 dss 3 (\B% + 2q 2 HlFh, 
Jo 

d 5 = d 4 - 3V2A / dss 3 Bl 
Jo 
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